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Abstract — Methods for primary user detection in cognitive 
radio may be severely impaired by noise uncertainty (NU) and 
the associated SNR wall phenomenon. The ability to avoid the 
SNR wall is proposed herein by detailed statistical modeling of 
the noise process when NU is present. A Gaussian model for the 
inverse noise standard deviation is proposed, and good agreement 
with the more common lognormal distribution is demonstrated 
for low to moderate noise uncertainty. Closed-form pdfs for a 
single noise sample and the energy of multiple noise samples 
are derived, allowing an optimal Neyman-Pearson detector to 
be employed when NU is present, thus avoiding the SNR wall 
effect. Initial measurements are presented that explore energy 
detection at low SNR in a practical system, showing that the noise 
distribution can be easily calibrated (learned) using a switch and 
matched load in the receiver. Useful detection performance down 
to -16 dB with energy detection is demonstrated, and it is found 
that noise uncertainty is not significant for an instrument-grade 
low-noise amplifier (LNA) for sub-minute acquisition times. 

I. Introduction 

Cognitive radio [1] is an interesting emerging paradigm 
for radio networks, where radios are able to sense and ex- 
ploit unused spectral resources, ideally improving spectrum 
utilization and allowing networks to operate in a more decen- 
tralized fashion. In the absence of cooperating primary users 
or beacons that indicate local spectrum usage, overlay-based 
cognitive radios must have sensing hardware and algorithms 
that are robust in the sense of providing very low missed 
detection rates at low SNR, thus impacting existing licensed 
users negligibly. 

Assuming an ideal noise model and given enough sensing 
time, simple methods like energy detection can theoretically 
discriminate the presence of a primary transmitter, even at 
very low SNR. In [2] the important effect of uncertainty in 
the noise distribution was identified and studied in detail, 
proving that when noise variance is confined to an interval but 
otherwise unknown, an SNR wall exists, below which useful 
detection performance cannot be guaranteed regardless of the 
observation time. 

The purpose of this paper is to study whether the SNR 
wall phenomenon can be eased by more detailed modeling 
of the noise uncertainty (NU) and to explore the impact of 
NU through direct measurement. It is argued in [2] that noise 
calibration to learn the noise model is not possible in cognitive 
radio since a primary may be present that corrupts any noise 
measurements. However, the dominant noise in RF systems 
is typically generated in the front-end amplifier and not the 



antenna. Thus, the noise level can be periodically measured by 
switching the receiver input to a matched load where primary 
signal is not present. Even then, the noise measurement may 
have error, but it is shown herein that with proper modeling 
of that error, the SNR wall can be reduced dramatically. 

The remainder of the paper is organized as follows. Sec- 
tion II provides some background on energy detection and re- 
views the SNR wall phenomenon. Section III defines the noise 
uncertainty model and derives closed-form noise pdfs in the 
presence of NU. Section IV provides numerical examples that 
illustrate how noise calibration can provide useful detection 
performance, whereas ignoring detection performance leads to 
an SNR wall. Section V presents initial measurements showing 
that energy detection at very low SNR is practically possible 
with the noise calibration technique. Finally, Section VI pro- 
vides some concluding remarks. 

II. Background 

This section briefly reviews concepts on optimal detection, 
energy detection, and noise uncertainty that are required for 
the remainder of the paper. 

The problem of primary detection in cognitive radio is 
usually treated using classical detection theory [3], where a 
decision must be made among two hypotheses: (Ho) only noise 
is present, or (Hi) signal plus noise is present. The received 
waveform x n under these two hypotheses is 

H 0 : x n = w n , n=l,2,...,N 
Hi : x n = w n + s n , n = 1, 2, . . . , N 

where w n and s n are the nth real noise and signal samples, 
respectively, and the detector must select Ho or Hi based only 
on observation of x n for n = 1,...,N. Given a decision 
rule, Pd is the probability of detection, or the probability 
that the detector correctly declares Hi, whereas Pf a is the the 
probability of false alarm, or the probability that the detector 
declares Hi when the true hypothesis is H 0 . 

When the pdfs of the received waveform x n under hy- 
potheses Ho and Hi are known, the Neyman-Pearson (N-P) 
detector provides optimal detection performance in the sense 
of providing maximum Pd for fixed Pf a . The N-P detector 
employs a likelihood ratio test (LRT), given by 
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where /ij(x) is the joint pdf of the observed samples for 
hypothesis H. For a selected threshold A, the detector declares 
Hi when £(x) > A, otherwise it declares H 0 . The threshold 
can be computed by fixing Pf a and inverting the cdf for the 
Ho (noise only) hypothesis. 

When noise and signal are both i.i.d. Gaussian, the energy 
of the signal, given by 



P 



N 
n=l 



(3) 



is a sufficient statistic. Detection based on p is known as 
energy detection, and the distribution of p is given by the 
Chi-Squared distribution, allowing the required LRT threshold 
and resulting detection performance to be computed in closed 
form. Given that the variance a 2 = Vax(w n ) is known, the 
energy detector can eventually provide near-perfect detection 
if N is made large enough, even at very low SNR. 

Unfortunately, a practical system will only have an estimate 
of Ctq, and this imperfect knowledge is referred to as noise 
uncertainty (NU). The NU concept was identified and studied 
in detail in [2], where noise variance is assumed to be confined 
to the interval [of 0 , a 2 { ] but otherwise unknown. In this case, 
worst-case detection performance for the N-P detector can be 
computed by assuming 



'hi' 



under H 



o- 



<7f 0 , under Hi, 



(4) 



thus providing the minimum separation of the H 0 and Hi pdfs. 
For a given noise interval, as the SNR is lowered a threshold is 
reached below which the worst-case energy detector exhibits 
-Pd < Pfa regardless of the number of samples. This complete 
detection failure is referred to as the SNR wall. 

III. Noise Uncertainty Modeling 

The main idea of this paper is to overcome the SNR 
wall phenomenon by more detailed modeling of the noise 
uncertainty. In this work, noise and signal are modeled as 
conditional Gaussian processes where a single real sample x n 
has the conditional distribution 



/(a 



2xp{-a 2 x 2 n /2}, 



(5) 



a = 1/ct, and a 2 is the variance. Note that the choice of 
using a rather than a as the modeled noise parameter in this 
work avoids having integration variables in the denominator, 
thus simplifying closed-form analysis. Given an i.i.d. process 
where a is fixed for a short time consisting of N samples, the 
marginal pdf of the vector x is 

If 00 f a 2 N 1 

/( x ) = QrfW J 0 f ^ aN CXP | T S j (6) 

where f(a) is the pdf of the unknown noise parameter a. Since 
the energy p — J2n=i x n * s a sufficient decision statistic here, 
we concentrate on this parameter. 



The distribution of p conditioned on a is given by the Chi- 
Squared distribution, or 

/(pI«) = 2^r( 2 jv/2) (a2p)A72 ~ 1 cx P{- a V2}, W 

and the marginal distribution f(p) therefore becomes 
1 



f(p) 



2 N / 2 T(N/2). 



Jo 



2 (a 2 p) N / 2 - 1 



cxp{— a p/2} da. 

(8) 

The idea of this paper is to choose a distribution for the inverse 
noise level f(a) that not only can be used to calibrate a 
practical system, but also has a simple form allowing (8) to 
be derived in closed form. 

The lognormal distribution is often proposed for modeling 
the variance of fading and noise processes, in which case f(a) 
is expressed as 

/lnO) = — 7=^= exp j-i(logCT-^ LN ) 2 /oL N l . (9) 

0-V27T(7 LN I 4 J 

where ^ln and ctln are the mean and standard deviation of 
log o\ Expressed in dB units ^ln = <Wb and ctln = SaaB, 
where 5 = log(10)/20. Letting a = 1 /a, (9) can be 
transformed to 

i'LJN(n) = ^==exp{ ^'logn • /M.\!" T\ S | 

' (10) 

which differs from (9) only in the sign of /^ln- A major draw- 
back of the lognormal distribution, however, is that closed- 
form analysis is often difficult. 

For small levels of noise uncertainty, we consider a much 
simpler model, where f(a) is assumed to be Gaussian which 
is fit to (10) using closed-form expressions for the mean and 
variance of (10) given by 

E {a} = cxpj-^LN + cr£ N /2}, (11) 

Std{a} = [cxp((Tln) - 1] cxp(-2^ LN + <7l N ), (12) 

where Std(-) denotes standard deviation. The pdf f(a) is then 
given by 



/(«) 



o, 



2na, 



expi-- 



1 (a - [i a ) 2 



a > 0, 



otherwise, 
(13) 

where the rescaling constant C a = crfc[— /i Q /(y / 2aJ)]/2 
results from the truncation of the left tail of the Gaussian at 
a = 0 and crfc(-) is the complementary error function. Note 
C a « 1 is omitted from later derivations, but it should be 
included if exact expressions are required. 

A. Single Sample: Marginal Distribution 

Assuming the Gaussian model for f(a), the marginal distri- 
bution of a single real sample is f(x) = J 0 °° f(a)f(x\a)da, 
or 

m ! 



2iro~ a 
1 

2na a 



ae -« 2 x 2 /2 e -( a -^) 2 /(2«l) da ^ (14) 



ae 



-[act 2 —ba+c] 



da, 



(15) 
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where 



respectively, followed by the substitution u = c\a 2 : 



x 2 1 

b = Ma/o£, 

c = £j{2a 2 a ). 



(16) 

(17) 
(18) 



The integral is of the form of the error function, which can 
be obtained by completing the square, resulting in 



/(*) 



47TC7 0 



e -cic 2 r^r 

h w— erfc(— v /cTc 2 )c2 

Ci V c l 



where 



ci = a, 



C2 = Ma/Va^ + 1), 



C3 = 



1^ 

2a2 



1 - 



1 



alx 2 + 1 



(19) 



(20) 
(21) 

(22) 



B. Multiple Samples: Energy Distribution 

For multiple independent samples, we will consider only 
the distribution of the energy p, which is a sufficient statistic 
whose distribution (8) becomes 



f(p) = 



1 



(23) 
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(24) 
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(26) 
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where 



a = ci 



2 + 2a£ 



c = /4/(2a*) 
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(30) 
(31) 

(32) 
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The integral I = h+ 12 can be evaluated by letting ii and 7 2 
be the contribution from a on the negative and positive axes, 



h = (-1) 7 
(-1) 



N - k I a N - k e- c ^ 2 da, 



(35) 



0 

N-k fC\c. 



2c^ k Jo 



t (JV-fc-l)/2 e -„ ^ (36) 



r(L fe ,ciC2)r(i fc ) 

where i fe = (TV + 1 - fc)/2, and 

r(a,x) = -^- [* e-H*- 1 dt (37) 
1 W Jo 

is the incomplete Gamma function. Similarly, 



I 



N-k-c x a 2 



2= a e 



da 



2cf 



-T(L k ). (38) 



Combining results in 

(39) 

C. Comparison of Gaussian and Lognormal 

It is instructive to consider in what situations the Gaus- 
sian assumption for f(a) provides a reasonable model. Fig- 
ure 1 plots f(a) side-by-side with f(x) (single sample) for 
^dB = 0 dB and different values of er d B € {0.5 dB, 1.0 dB, 
2.0 dB}. A log scale is used for f(x) to highlight the small 
differences in the distribution tails. For small and moderate 
levels of noise uncertainty, the Gaussian approximation for 
f(a) is very close to the lognormal model. Also, for low noise 
uncertainty, the small mismatch in f(a) results in negligible 
error in the marginal density f(x). For larger noise uncertainty, 
significant differences in the two models can be seen. 

IV. Detection With NU 

In this section we demonstrate with a simple example how 
having a model of the noise uncertainty can increase detection 
performance and remove the SNR wall. In this example, 
a is considered to be an unknown parameter following a 
lognormal distribution with h^b —0 dB and CTdB=l dB, which 
is subsequently fit using a Gaussian distribution. Signal and 
noise variance are assumed to be equal (SNR=0 dB) and 
N = 20 samples are used for detection. 

First, a worst-case analysis like that presented in [2] is 
considered. Here, only bounds are set on the noise level, and 
the structure of the noise variation is ignored. It is assumed 
that the worst-case values for a are \i a ± 1.5a a , which is 
conservative since the a will sometimes fall outside of these 
bounds. Figure 2 shows the Chi-Squared pdfs for the worst 
case assumption (4), indicating that detection is not possible 
since the Ho curve is actually to the right of the Hi curve. 

Next, the structure of the noise error in (39) is taken into ac- 
count, producing the pdfs in Figure 3 and indicating sufficient 
separation for useful detection. Figure 4 shows the detection 
performance from the worst case analysis and the case that 
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Fig. 1. A comparison of the lognormal distribution on a and a Gaussian approximation of the same distribution for different levels of noise uncertainty: 
°"ds G {0-5 dB, 1.0 dB, 2.0 dB}. The noise uncertainty pdf f(ct) is plotted on the left, and the corresponding marginal single sample pdf f(x) is plotted 
on the right for each level of noise uncertainty. MC gives the results of Monte-Carlo simulations of the exact distribution with lognormal NU, compared with 
the Gaussian NU approximation (Gauss), and no noise uncertainty (No NU). 



exploits the noise error pdf. In the worst-case analysis, the 
SNR wall has clearly been crossed, since Pd < Pf&- On the 
other hand, exploiting the known statistics of the noise error 
allows useful detection even when the exact noise level is 
uncertain. 

V. Noise Calibration Measurement 

In this section we present the results of an experiment that 
tests the possibility of energy detection at low SNR using 
practical hardware. As indicated in [2], noise calibration can be 
difficult in traditional wireless receivers where it is unknown if 



the primary is present or not. However, since most of the noise 
in a true receiver comes from the front-end low-noise-amplifier 
(LNA), the simple architecture depicted in Figure 5(a) can be 
used for noise calibration. To learn the noise distribution, the 
cognitive radio node periodically switches the receive channel 
away from the antenna to the matched termination to sample 
and learn the noise distribution. 

This idea was tested using the experimental setup shown 
schematically in Figure 5(b). The setup employs a custom 
multiple-input multiple-output (MIMO) channel sounder that 
is basically equivalent to that presented in [4], with the 
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Fig. 2. Distribution of f(p) for noise (Ho) and signal plus noise (Hi) a t> u u-r t c j . n u u i » t-i- 1 1 r> f 
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Fig. 3. Distribution of f(p) for noise (Ho) and signal plus noise (Hi) using 
the proposed Gaussian NU model 



exception that custom FPGA-based data acquisition is used 
in the present system. 

The transmit (TX) node simulates the primary user, where a 
baseband Gaussian signal with a flat W = 20 MHz bandwidth 
is generated in 100 /is frames with the arbitrary waveform 
generator (AWG), up-converted to 2.55 GHz, power amplified 
to 23 dBm, and fed to either the active transmit channel (TX1) 
or a matched load (TXO). The channel is a simple direct cable 
connection from the transmitter to receiver, where different 
fixed attenuators are inserted giving loss L and producing 
different SNR levels at the receiver. 

The receive (RX) node simulates the cognitive radio that 
employs a switch to feed its single receive chain either from 
the channel (RX1) or from a matched load (RXO). The receive 
chain consists of a 40 dB wideband LNA, down-conversion 
to a 50 MHz IF, followed by FPGA-based f s = 200 MS/s 
data-acquisition. For this experiment, the raw IF samples are 
stored, passed to a PC, down-converted to complex baseband, 
and filtered (20 MHz bandwidth) using MATLAB before 
performing energy detection. 

A total of M data records are acquired during each mea- 
surement, where the mth record is depicted in Figure 5(c). 
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Fig. 5. Measurement setup for experimental study: (a) envisioned cognitive 
radio employing noise calibration, (b) channel sounder based acquisition 
system for experiment, (c) acquisition frame structure 



Four phases are used in each record to probe all four switch 
combinations. Within a single phase, the channel is acquired 
for T — 100 /is followed by a delay of To, where To = 0 can 
be used for back-to-back acquisition. During post-processing, 
only N s samples within each acquisition window are used, 
thus spanning time T s = N s / f s in order to simulate different 
integration windows in a cognitive radio energy detector. We 
will denote the nth filtered complex-baseband sample, of the 
kth phase, in the mth record as x m ,k, n - The four phases are 
denoted symbolically as k e{TXl RX1, TXO RX1, TX1 RXO, 
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Fig. 6. Empirical pdfs (solid lines) for the four measurement phases 
with parameters: T D =0, L=110 dB (p = -6 dB), N s = 100, and 
M = 800 realizations. Fitted stationary Chi-Squared distributions are also 
shown (dashed lines). 



TXO RXO}. 

The energy in each record is computed using a variable 
window size N s according to simple integration, or 

Pm,k = ^ l^m.fc.nl 2 , (40) 
n=iVi+l 

where N\ = 50 samples are always skipped at the beginning 
of each frame to avoid artifacts from the switching operations. 
Empirical distributions of the noise and signal plus noise 
energy are finally computed with a histogram using the M 
energy snapshots in (40). 

Empirical signal and noise pdfs are also compared with ideal 
Chi-Squared pdfs with N — 2WN s /f s degrees of freedom 
with sample variance estimated according to 



Fig. 7. Probability of detection (Pd) versus probability of false alarm (Pf a ) 
for To = 0, L = 110 dB (p = -6 dB), and N s e{100,800,6400}. 
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i 



M 



Pn 



(41) 



The SNR (p) in dB is estimated at an attenuation level of 
L =80 dB (high SNR) using 



p(L = 80 dB) = 10 log! 



'TX1RX1 



(42) 



, CT TX1RX0 - 

and SNR at lower attenuation levels is computed using p(L) = 
p(L = 80 dB) + 80 — L. Note that for convenience, pdfs 
are normalized with respect to the expected noise energy 

Figure 6 plots empirical noise/signal pdfs for the four phases 
compared with Chi-Squared pdfs for To = 0 (back-to-back ac- 
quisition), L — 110 dB (p — —6 dB), and N s = 100 samples 
(N = 40). First, the simple Chi-Squared distribution (no noise 
uncertainty) provides a good fit to the empirical pdfs for all 
cases, suggesting that the noise parameter a is fairly constant 
over the total acquisition time of ATM = 320 ms. Also, there 
is no apparent difference in the energy pdfs for the three noise- 
only (Ho) phases, indicating that having the RX connected to 
a matched load (TX1 RXO and TXO RXO) is equivalent to 
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Fig. 8. Empirical pdfs (solid lines) for the four measurement phases with 
parameters: T D =24 ms, L=110 dB (p = -6 dB), N s = 100, and M = 
600 realizations. Fitted stationary Chi-Squared distributions are also shown 
(dashed lines). 



measuring the channel with the transmitter not present (TXO 
RX1). Finally, good separation of the pdfs for the Ho phase 
and Hi phase (TX1 RX1) is demonstrated, indicating that 
useful detection with the direct noise measurement is possible. 

Figure 7 plots probability of detection versus probability 
of false alarm Pf a for the same case, computed directly 
from the empirical pdfs, where phase TX1 RXO is used to 
estimate hypothesis Ho. The result shows that when the noise 
distribution is measured, near perfect detection is possible if 
the sample size is made large enough. 

Figure 8 plots the empirical pdfs of the four phases for 
the same case, but with a longer delay between acquisitions 
To = 24 ms. Although we expected that increased noise 
uncertainty would result and spread the pdfs for the longer 
acquisition time of A{T D + T)M = 58 s, they still exhibit 
an excellent fit to the simple Chi-Squared distribution with 
no noise uncertainty. Possible reasons that the noise process 
is so stable is that an instrument-grade LNA is used and the 
temperature in the channel sounder was likely very constant. 
In the future, we intend to study inexpensive consumer-grade 
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Fig. 9. Empirical pdfs (solid lines) for the four measurement phases with 
parameters: T D =24 ms, L = 120 dB (p = -16 dB), N a = 6400, and 
M = 600 realizations. Fitted stationary Chi-Squared distributions are also 
shown (dashed lines). 




Fig. 10. Probability of detection (Pd) versus probability of false alarm (Pf a ) 
for T D = 0, L = 120 dB (p = -16 dB), and N a e{ 100,800,6400} samples 



LNAs under varying environmental conditions, where noise 
uncertainty is more likely to occur and noise calibration is 
more challenging. 

Figure 9 shows pdfs for the case of lower SNR L = 120 dB 
(/> = -16 dB), T D = 24 ms, and N s = 6400 samples. 
Even though the separation of the pdfs is poorer in this case, 
the detection performance plotted in Figure 10 indicates that 
the SNR wall effect is avoided, and useful detection is still 
possible with a long enough integration time. 

In the future, the system will be modified to allow longer 
integration times where the effects of noise uncertainty are 



more likely to be observed. However, these results suggest 
that detection at SNRs of -16 dB or lower should be possible 
using simple noise calibration. 

VI. Conclusion 

Although noise uncertainty can be a severe impairment to 
robust detection in cognitive radio at low SNR, this paper 
has shown that the SNR wall effect can be overcome by 
proper specification of the noise variation. A closed-form pdf 
that assumes a Gaussian distribution for the inverse noise 
standard deviation was derived and it was shown that the 
model provides a good fit to the more commonly assumed 
lognormal pdf for low to moderate noise uncertainty. A 
simulation example was presented, confirming that by properly 
modeling the noise uncertainty, the SNR wall phenomenon can 
be avoided, providing useful energy detection performance at 
very low SNR. 

Initial experimental measurements were also presented that 
explore energy detection performance in a true receiver using 
practical hardware. Detection performance based on empir- 
ically measured pdfs indicated that useful detection down 
to at least -16 dB is possible with energy detection using 
sufficient integration time. Measured noise distributions over 
short (0.3 s) and moderate (58 s) acquisition times showed 
negligible deviation from a Chi-Squared distribution, suggest- 
ing that the noise level in our system is very stable and that 
detailed modeling of the noise uncertainty is unnecessary for 
sub-minute integration times. 

Since only an expensive instrument-grade LNA was consid- 
ered in this work, future work will explore noise variation in 
low-cost commercial-grade amplifiers that may exhibit noise 
statistics that are much less stable. Additionally, due to existing 
limitations of the system acquisition firmware, we could not 
explore the performance of very long acquisition times, which 
will also be the subject of future investigations. 
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